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The phenomenological two-level atom is re-analysed using the methods of effective field the- 
ory. By presenting the Dicke-Jaynes-Cummings model in real space, an exact diagonalization is 
accomplished going beyond the rotating wave approximation. The meaning of the symmetries and 
conserved quantities in the theory is discussed and the model is related to non-equilibrium field 
theory. The structure of the solution raises a question about the rotating wave approximation in 
quantum mechanics. 



It was first realized by Jaynes and Cummings that a semi-classical model of a two-level atom could reproduce the 
essential features of the quantum theoretical problem . Recently there has been a ressurgence of interest in effective 
field theory in several works particularly in connection with quantum electrodynamics |^-^ ■ Effective theories are often 
the most efficient way of performing calculations in field theory; they naturally separate conceptually independent 
issues. An alternative to an effective theory is to perform ab initio calculations, first principle computations in 
^ : quantum electrodynamics, but this is orders of magnitude harder and involves computing effects in which one is not 
interested. For instance, if one is not interested in the reason for the energy spacing of the bound state levels in the 
atomic system, it is not neccesary to derive it: rather it can be inserted either from experiment of other calculations. 
' This note reexamines an extremely well known and widely used model for the resonant interaction between a two 
I level atom and radiation. The two level system has a broad repertoire of applications in physics, from spin models to 
00 ■ the micromaser It is related to a class of Dicke models and, in the so-called rotating wave approximation, it 
0^ ' becomes the Jaynes- Cummings model ^ which may be solved exactly. A Hamiltonian analysis of symmetries in this 
, Jaynes-Cummings model is given in ref. In reference |^ a variation on Schwinger's closed time path | |To[ | method 
^ ^' of analysing non-equilibrium systems was presented in which field theories with spacetime dependent perturbations 
^ I , ,, could be formulated analogously to gauge theories. The same technique can be applied to the two-level atom to solve 
(~| ■ the full model and eliminate the need for the so-called rotating wave approximation. In doing so, one obtains extra 
Ci ' physical insight into the system as well as a result which is valid over a wider range of parameters. It is also possible 
^ • to verify a previous claim in ref. Q about the relationship between non-equilibrium effective field theory and Rabi 
oscillations of the two-level atom. 

Consider the phenomenological two-level system described by the action 



S= dVt 
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where A,B — 1,2 characterizes the two levels, ihDt = ihdt + i^{t) in matrix notation, and F — Tj^g is an off- 
diagonal anti-symmtrical matrix. At frequencies which are small compared to the light-size of the atom, an atoms 
may be considered electrically neutral. The distribution of charge within the atoms is not required here. In this 
approximation the leading interaction is a resonant dipole transition. The connection F^^ plays an analogous role to 
the electromagnetic vector potential in electrodynamics, but it possesses no dynamics of its own. Rather is works as a 
constraint variable, or auxiliary Lagrange mulitplier field. There is no electromagnetic vector potential in the action 
since the field is electrically neutral in this formulation. Tab refers not to the U{1) phase symmetry but to the two 
level symmetry. Variation of the action with respect to T{t) provides us with the conserved current. 

-irA^s-rsM (2) 
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which represents the amplitude for stimulated transition between the levels. The current generated by this connection 
is conserved only on average, since we are not taking into account any backreaction. The conservation law corresponds 
merely to 

dt oc sin(2 1 Xit)) (3) 
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where X{t) will be defined later. The potential VAsit) is time dependent and comprises the effect of the level 
splitting as well as a perturbation mediated by the radiation field. A 'connection' = — is introduced since 
the diagonalization procedure requires a time-dependent unitary transformation and thus general covariance demands 
that this will transform in a different basis. The physics of the model depends on the initial value of this 'connection' 
and this is the key to the trivial solubility of the Jaynes-Cummings model. 
In matrix form we may write the action for the matter fields 

dVt rAOABi^B (4) 

where 



o = 



2m ' ^ 2 



Jit) - -V2-fDt 



(5) 



The level potentials may be regarded as constants in the effective theory. They are given by Vi = Ei and V2 = E2 — Mlii 
where hQ.ii is the interaction energy imparted by the photon during the transition i.e. the continuous radiation pressure 
on the atom. In the effective theory we must add this by hand since we have separated the levels into independent 
fields which are electrically neutral; it would follow automatically in a complete microscopic theory. The quantum 
content of this model is now that this recoil energy is a quantized unit of hil, the energy of a photon at the frequency 
of the source. Also the amplitude of the source J would be quantized and proportional to the number of photons 
on the field. If one switches off the source (which models the photon's electric field) this radiation energy does not 
automatically go to zero, so this form is applicable mainly to continuous operation (stimulation). The origin of the 
recoil is clear however: it is the electromagnetic force's interaction with the electron, transmitted to the nucleus by 
binding forces. What we are approximating is clearly a J**^^ term for the electron, with neutralizing background 
charge. 

It is now desirable to perform a unitary transformation on the action ^jj Utp, O — > UOU~^ which diagonalizes 
the operator O. Clearly the connection Tab will transform under this procedure by 

?'?) 

T^T+-{U{dtU-')-{dtU)U-') (6) 

since a time-dependent transformation is required to effect the diagonalization. For notational simplicity we define 
L = — — Df, so that the secular equation for the action is: 

{L-Ei- A)(L -E2 + nn-x)- (j2 + rj^) = o. (7) 

Note that since J dt J = there are no operator difficulties with this equation. The eigenvalues are thus 



X±= L- Ei2 + hQ±\l^ {E21 - + J2 + TI2 (8) 



= L- Ei2 + 'h^l± ^Jn^Cb'^ + J2 + TI2 (9) 

= L-Ei2+h9.± huR, (10) 

where E12 = ^(Ei + E2) and E21 = {E2 — Ei). For notational simplicity we define oj and uir. One may now confirm 
this procedure by looking for the eigenvectors and constructing as the matrix of these eigenvectors. This may be 
written in the form 

^_i^/cosJ -sme\ 

\ smO cos 6 I ^ ' 



where 



cos 6*= — -j^=^^^=^^^=^= (12) 



^n'(u; + a;ii)2 + j2 + r22 



sin^= , V-^'+I^i2 (13) 



^n2(u; + a;,i)2 + J2 + r22 
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The change in the connection r{t) is thus off-diagonal and anti-symmetric as required by the gauge symmetry con- 
servation law: 



udfU- 



dtO 

-dtO 



The time derivative of 9{t) may be written in one of two forms which must agree 

{dtO) = = ^• 

— sin cos 

This provides a consistency condition which may be verified and leads to the proof of the identities 

ojRdtOJR = JdtJ + TdtT 
Vj2 + r2(9t + A)v/j2+r2 + (a) + WR){dt + A)(a) + wr) = 

for arbitrary J{t) and T{t), where 

1 dt{{co + LUR)' + j^ + r') 



(14) 



(15) 



A 



2 {Q + ur)^ + J2 + r2 



(16) 



(17) 



These relations are suggestive of a conformal nature to the transformation and, with a little manipulation using the 
identities, one evaluates 



{JdtJ + TdtT) 
OjRy/J-'+r' 



1 



{W + U}r){Q + 2llIr) 

(a) + a;fl)2 + j2 + r2 



(18) 



This quantity vanishes when -|- is constant with respect to time. Owing to the various identities, the result 
presented here can be expressed in many equivalent forms. In particular, it is zero when a) = 0. The equations of 
motion for the transformed fields arc now 



£^12 + hu)R idtO 
—idtO L — E12 — hujR 



= 0. 



(19) 



In this basis, the centre of mass motion of the neutral atoms factorizes from the wavefunction, since a neutral atom 
in an elect romagntic field is free on average. The two equations in the matrix above may therefore be unravelled by 
introducing a 'gauge transformation' or 'integrating factor' 



where the free wavefunction in n ■ 



3 dimensions is 



pi(k-x— cut 



(20) 



(21) 



for the centre of mass motion is a general linear combination of plane waves satisfying the dispersion relation for 
centre of mass motion 



2 1,2 



2m 



+ h{Q-ij)-E 
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0. 



(22) 



The latter is enforced by the delta function. This curious mixture of continuous (lu) and discontinuous (fi) belies the 
effective nature of the model and the fact that its validity is only for a continuous operation (an eternally sinusoidal 
radiation source which never starts or stops). The relevance of the model is thus limited by this. Substituting this 
form, we identify X{t) as the integrating factor for the uncoupled differential equations. The complete solution is 
therefore 



(23) 



Notice that this result is an exact solution in the sense of being in a closed form. In the language of a gauge theory this 
result is gauge dependent. This is because our original theory was not invariant under time dependent transformations. 
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The covariant procedure we have apphed is simply a method to transform the equations into an appealing form; it 
does not imply invariance of the results under a wide class of sources. 

That this system undergoes transitions in time may be seen by constructing wavefunctions which satisfy the bound- 
ary conditions where the probability of being in one definite state of the system is zero at t — 0. To this end we write 
^1 = + V'-) E^nd — ~ V-"-)- In order to proceed beyond this point it becomes necessary to specify the 

initial value of ri2. This choice carries with it physical consequences; the model is not invariant under this choice. 
The obvious first choice is to set this to zero. This would correspond to not making the rotating wave approximation 
in the usual two level atom, with a cosine perturbation. Focusing on the state \l/o which was unoccupied at t = for 

ri2 = 0, 




^ Jofi sin(m') 



2hujR 



^j{x). (24) 



We are interested in the period and amplitude of this quantity, whose squared norm may be interpreted as the 
probability of finding the system in the prepared state, given that it was not there a,t t — Q. Although the integral 
is then difficult to perform exactly, it is possible to express it in terms of Jacobian elliptic integrals, logarithms and 
trig functions. Nevertheless it is clear that Co = \{E2i/h — 51) is the decisive variable. When huj <^ Jq is small, 
the first term is Jq cos{ilt) and the second term is small. This is resonance, although the form of the solution is 
perhaps unexpected. The form of the wavefunction guarantees a normalized result which is regular ai & — and 

one has ~ sin ^ dt' ^ cos{flt')j , which may be compared to the standard result of the Jaynes-Cummings model 

^0 ~ sm{Jot/h). In the quantum case the amplitude of the radiation source Jq is quantized as an integral number Nq 
of photons of frequency Q. Here we see modulation of the rate of oscillation by the photon frequency (or equivalently 
the level spacing). In a typical system the photon frequency is several ten orders of magnitude larger than the coupling 
strength Jq <gC Ufl ~ E12 and thus there is an extremely rapid modulation of the wavefunction. This results in an 
almost chaotic collapse-revival behaviour with no discernable pattern, far from the calm sinusiodal Rabi oscillations 
of the Jaynes-Cummings model. If huj ~ Jq the second term is of order unity and then, defining the normalized 
resonant amplitude 

, ^° (25) 



one has 



^„ ^ sin I ^ £; {Qt, A)-A [ d{nt) ''"^"^^ - I V^(x). (26) 



A 



1-A^ sin2(m) 



The Jacobian elliptical integral E{a,P) is a doubly periodic function so one could expect qualitatively different 
behaviour away from resonance. On the other hand, far from resonance huj 3> Jq, the leading term of the connection 
becomes ^"0 ^ sin {Cjt) 1/1(2;) ~ sin (Qt) 'ip{x), and the effect of the level spacing is washed out. 

One can also consider other values for the connection. Comparing to the off diagonal sources ^^D^, predicted 
on the basis of unitarity in effective non-equilibrium field theory |9), one obtains an indicatation that, if the initial 
connection is in phase with the time derivative of the perturbation, then one can effectively 'resum' the decay processes 
using the connection. This is a backreaction effect of the time dependent perturbation, or a renormalization in the 
language of ref. Q. If one chooses = Jo sin(rit), this has the effect of making the off-diagonal terms in the action 
not merely cosines but complex a conjugate pair Jq exp(±ir2<). This corresponds to the result one obtains from 
making the rotating wave approximation near resonance. This initial configuation is extremely special. With this 
choice, one has exactly 



*o = sin / dt' 



^(x). (27) 



The stability of the solution is noteworthy, and the diagonalizing transformation is rendered trivial. The connection 
dtO is now zero under the diagonalizing transformation. Thus the above result is exact and it is the standard result 
of the approximated Jaynes-Cummings model. This indicates that the validity of the Jaynes cummings model does 
not depend directly on its approximation, but rather on the implicit choice of a connection. 

The Jaynes-Cummings model is at the heart of many laser phenomena. Korenman attempted to create a detailed 
microscopic non-equilibrium field theory of the laser from the Lagrangian above; a recent letter ||l^ proposes to 



4 



do a similar thing for semiconductor lasers. In the first of these Korenman claims to present a microscopic theory, 
but the dipole interaction he considers is an effective interaction, as are the variables for the pump part of the laser, 
so it is clearly not an ab initio approach from the level of electron interactions. Nevertheless, the structural elements, 
including exchange mechanisms are present. It is interesting to see how the model presented here is related to the 
non-equilibrium field theory approach presented in ref. Q , which is directly related to the method used by Korenman. 
There is a direct analogy between the two level system and the forward and reverse time directions of the Schwinger 
closed time path method, dictated by unitarity. Moreover the generalized model presented in ref. contains the 
present model in the non-relativistic limit. 

The form of the action in eqn (1) appears arbitrary but it may be understood in the context of ref. ||^ a conformal 
perturbation. The action is perturbed by a time-dependent source which one hopes would lead to a stable theory (the 
form of the action remaining the same over time). This suggests conformal covariance. The conformal link can be 
made by writing the above model as the limit of a pseudo-rclativistic theory since the conformal group is an extension 
of the Lorentz group. This also makes contact with ref. The consistency of such an approach has been verfied in 
ref. [T^ . Beginning with the Lorentz covariant action 

J dVtf^^{df,c^A){d^M + lrnA(l)\ + (j>AJAB(l}BY (28) 

we consider a conformal rescaling by letting 5^,^ il^ 'g^^. The action is not invariant under this rescaling: if it were, 
there would be no need for the connection F, or indeed this paper. The volume element scales as the square root of 
the determinant of the metric, i.e. ^ — > fi'^y^ in 3-1- 1 dimensions, but we shall keep this separate for now. Since the 
issue is not invariance but equivalence, this will not play a crucial role. The first term in the braces contains one inverse 
power of the metric, the second none and the third two. Choosing fl^ = Jab, the off-diagonal, symmetric matrix 
with non-zero elements J(i), one can absorb the time dependent interaction by performing a generalized rescaling. 
Rescaling the fields by ^ il4>, the action takes the form 

S = J dVl \^-{D;4>A)iD^l>A) + I^'a^a - i1b1'{d;4>)^ , (29) 

where — O^Sab + Iab^ which is obtained by moving the scale factors through the derivatives, and 

IdtJAB fdtVL\ 

lAB- (30) 
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The familar form of the conformal correction d^il/VL is replaced in eqn. (1) simply by T12 — dt^, which makes the 
initial value of ^ dtJ clear: it is the connection required for the derivatives to commute with a conformal rescaling 
brought about by the perturbation. 

The non-relativistic limit of the transformed action with a time-only dependent U,{t) leads to eqn. (1), up to a 
Jacobian. Thus although these actions are not identical, they are related by an overall factor which behaves as though 
to view the system through a distorting glass. A similar non-equilibrium situation may be found in ref. |p^ . These 
two theories give essentially the same results because they have the same structural elements. Indeed, it might be 
possible to show that they form an equivalence class in the formal sense, though this has not been shown here. It 
is this feature of effective field theories which makes them robust and usable. The result of this analysis is perhaps 
surprising though: the result which one obtains by making the rotating wave approximation ~ sin(r2t) in a 
quantum mechanical formalism (the JC model) is a stable, steady state effective field theory, but the same model 
solved fully without approximation is not stable to rescalings and leads to extremely complicated behaviour. This 
is due as much to the method of analysis in time dependent quantum mechanics as it is to do with the rotating wave 
approximation. Here we have an example where such a model can be solved exactly and does not appear to be give 
any qualitative pattern of behaviour at experimental values typical for the micromaser Q . What is interesting is why 
the specific experimental situation should be well described by the very special Jaynes-Cummings model rather than 
one of the many less stable theories with a different initial ri2. This certainly warrants further investigation, perhaps 
in the context of conformal field theory. 
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